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An optimal estimator for the CMB-LSS angular power 
spectrum and its application to WMAP and NVSS data 
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ABSTRACT 

We use a Quadratic Maximum Likelihood (QML) method to estimate the angular 
power spectrum of the cross-correlation between cosmic microwave background and 
large scale structure maps as well as their individual auto-spectra. We describe our 
implementation of this method and demonstrate its accuracy on simulated maps. We 
apply this optimal estimator to WMAP 7-year and NRAO VLA Sky Survey (NVSS) 
data and explore the robustness of the angular power spectrum estimates obtained 
by the QML method. With the correction of the declination systematics in NVSS, 
we can safely use most of the information contained in this survey. We then make 
use of the angular power spectrum estimates obtained by the QML method to derive 
constraints on the dark energy critical density in a flat ACDM model by different 
likelihood prescriptions. When using just the cross-correlation between WMAP 7 year 
and NVSS maps with 1.8° resolution, the best- fit model has a cosmological constant 
of approximatively 70% of the total energy density, disfavouring an Einstein-de Sitter 
Universe at more than 2 a CL (confidence level). 

Key words: cosmic microwave background - large scale structure - methods: numer- 
ical - methods: statistical - cosmology: observations 



1 INTRODUCTION 

Understanding of the nature of dark energy is one of the 
outstanding questions in observational cosmology. Since the 
discovery of the present acceleration of the Universe by 
the measurement of the luminosity distance of distant type 
la su pernova (SN la) (|Riess et al.l 119981: IPerlmutter et al.l 



19991'). several obser v ations (e.g., Tegmark et al.l |2004| : 



Eisenstein et al.l uOOSi : iLarson et al.l uOllf ) have converged 



to a cosmological concordance model in which an unknown 
component having a negative pressure density - 'dark en- 
ergy' - that contributes ~ 2/3 of the total energy budget of 
the Universe. At present the precise nature of dark energy, 
parameterised by its equation of state, can only be weakly 
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constrained using a range of cosmological tests, but indica- 
tions are that its behaviour is close to that expected from a 
cosmological constant. 



A key strategy in determining the nature of dark energy 
is to combine as many different observations as possible, in- 
cluding the luminosity distances for Type la supernovae, the 
baryonic acoustic oscillation (BAG) scale observed in galaxy 
surveys, anisotropies of the cosmic microwave background 
(CMB) and weak lensing surveys. Cross-correlations among 
the above observations also contain precious cosmological 
information about dark energy. Ambitious space projects 
have been proposed to address the dark energy questio n 
with this strategy, including EUCLID (jLaureiis et al.ll2009l ). 
which will focus on BAG and weak lensing, and WFIRST 
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Q an infrared satellite with a focus yet to be specified. 
In the meantime, ground based programs such as DES Q 
PanSTARRS0, LSStEI will also improve the current under- 
standing of structure formation and provide excellent galaxy 
surveys to cro ss-correlate with the CMB anisotropy maps 
from Planck (|Planck Collaboration! 120061 ). 

One of the key indicators of the presence of dark 
energy are CMB fluctuations created by the l ate In tegrated 
Sachs- Wolfe (ISW) effect (|Sachs fc Wolfd 1 19671 ). When 
the Universe is not completely matter dominated, CMB 
anisotropics are created at l ate times and these contribut e 
most at large angular scales (jKofman fc Starobinskvlll985h . 
Since the low multipoles of the CMB angular power 
spectrum are mostly affected by cosmic variance, an 
extraction of the ISW part solely from CMB data is rather 
difficult, but it is feasible when CM B is cross-correlated 
with large scale structure (LSS) (jCrittenden fc Turokl 
[l995|). Several positive detections of the ISW-LSS cross- 
correlation have been performed since the release of the 
WMAP first year data by using different tracers of LSS 
and s t atistical estimators ( e .g. Afshordi. Loh fc StraussI 
p004): iBoughn fc CrittendenI (|2004|) : iFosalba et al.l (|2003|) 



iNolta at al. I (|2004): IVielva. Martinez-Gonzalez fc Tucci 



200e): IPictrobo n. Balbi and Marinuccil (|2006|): IHo et a" 
200g ): iGiannantonio et al.1 (|2008l ): iFrancis fc Peacock 



2010 ). lDupe et al.l ( 201ll ) and references therein) 



One of the purposes of this paper is to develop tools 
to estimate the angular power spectrum (APS) of the cross- 
correlation between CMB and LSS by a quadratic maximum 
likelihood (QML) method. The QML method in this con- 
text has a number of advantages: foremost, given the low 
signal-to-noise expected for the ISW-LSS cross-correlation, 
it is essential to use a minimum variance method, such 
as QML, to estimate the cross power spectrum. In addi- 
tion, based in pixel space, the QML method is ideal for 
accounting for the incomplete sky coverage and masks of 
the surveys. Finally, while the QML method is expensive 
computationally, the fact that the ISW signal is primar- 
ily at low multipoles means that it is tractable to con- 
strain it on maps using only a modest resolution. The 
QML method has also found application in the estima- 
tion of the power spectru m of the CMB intensity and 
polarization (Tcgmark 19971; iTegmark fc de Oliveira-Costal 
I2OOII ) and has been re cently applied to the latest re- 
lease s of WMAP data llCruppuso et al.l l2009l : iPaci et al.1 
l2010l : lGruppuso et al.ll201ll ). A QML estimator was already 
used to measure th e CMB -LSS cross-correlation only by 
IPadmanabhan et al.l (|2005j), however our implementation is 
different in few important aspects: the inversion of the matri- 
ces is implemented here using the single value decomposition 
(see also Section 3.1) and all the three spectra - TT, TG, GG 
- are computed for all the multipoles in the range of interest. 

Another purpose of this work is to apply our method- 
ology to availabl e public CMB and LSS data, namely 
WMAP 7 year djarosik et al.1 I2OIII) and N RAO VLA 
Sky Survey (NVSS) data (|Condon et allll998l ). NVSS has 



been one of the most widely used surveys in the con- 
text of ISW studies because the radio galaxies it sur- 
veys are at high redshifts and it covers a large sky frac- 
tion of the sky; however, contradicting claims about the 
evidence of its non- vanishing correlation with CMB ex- 
ist in the literature ([Pietrobon. Balbi and Marinuccil I2OO6I : 



ISawangwit et al.ll2010l ) (see also iDupe et al.l (2011) for an 
exhaustive compilation of existing results). It is therefore 
important to apply an optimal methodology to address and 
quantify the evidence of cross-correlation between the most 
recent large scale CMB measurement and one of the largest 
LSS survey available. 

Our paper is organized as follows: in Section [2] we de- 
scribe the QML method and give technical details of our 
implementation of it; Section [3] discusses our tests of the 
implementation in simulated maps. In Section!?] we report 
the APS estimates obtained from WMAP 7 year and NVSS 
data, and then we use these estimates of the cross-correlation 
in Section[S]to derive constraints on the present critical den- 
sity due to the cosmological constant. Finally in Section [6] 
we draw our conclusions. 



2 METHODOLOGY 

2.1 The QML approach 

The quadratic maximum likelihood method for power spec- 
trum est i mate of CMB anisotropies was introduced by 
TegmarkI (119971) and later extende d to polarization by 



Tegmark fc de Oliveira-Costal (|200ll ). Previously, a QML 



was employed to me asure t he cross-correlation between 
CMB and LSS only bv IPadmanabhan ct al. (200i) (see also 
IHo et al.l (|2008l )). The code in IPadmanabhan et al.1 (|2005l ) 
estimated only the cross-correlation power spectrum only, 
with a fast and approximated algorithm to invert matri- 
ces and used the approximation of a block diagonal covari- 
ance matrix. In what follows we shall describe the QML 
method for the whole CMB-LSS data and our implementa- 
tion wh ich does not depend on the simplifying assumptions 
used in IPadmanabhan et al.l (|2005h . 

Given a CMB map in temperature and a galaxy survey 
x= (T, G) (vector in pixel space), the QML provides an 
estimator of the angular power spectrum Cf - with X being 
one of TT, TG, GG. This estimator is given by 



C 



E(^- 



^XX' 

)U' 



X E^/ X - 



-tr(NEf) 



where the F, 



XX' 



w 



is the Fisher matrix defined as 



r,XX' 1, 
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and the E matrix is given by 
1 dC ^_i 
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2 dC'^ 



(1) 



(2) 
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^ http://wfirst.gsfc.nasa.gov 

^ http://www.darkenGrgysurvey.org/ 

^ http://pan-starrs.ifa.hawaii.edu/public/ 

* http://www.lsst.org/lsst/ 



C = S{Cf) + N being the total global covariance matrix 
including the signal S and noise N contributions. Cf is 
called the fiducial theoretical power spectrum and also is 
used to create the simulated maps used to test the method 
in Sec. |3] 
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Although an initial assumption is needed for this fidu- 
cial power spectrum, the QML method provides unbiased 
estimates of the power spectrum of the map regardless of 
this initial guess 



(Cf) = Cf 



(4) 



Here the average is taken over the ensemble of realizations 
based on the input spectrum Cf . (See Sec. |3]for more de- 
tails.) The assumed fiducial power spectrum can impact the 
error estimates, but in practice we start near enough to the 
true result to be able to neglect this effect. The QML method 
is also optimal, since it can provide the smallest error bars 
allowed by the Fisher-Cramer-Rao inequality. 



X'\ 



{AC^ACP > = (F 



-i^x X' 



where 



{Cf), 



(5) 



(6) 



and the averages, as above, are over an ensemble of realiza- 
tions. 

Our implementation of the QML method is fully par- 
allelized (MPI) and written in Fortran 90. The inversion of 
the covariance matrix C scales as 0{N'^^^). The number of 
operations is roughly driven, once the inversion of the total 
covariance matrix is done, by the matrix-matrix multiplica- 
tions to build the operators E^ in Eq. (|3} and by calculat- 
ing the Fisher matrix F^i^ given in Eq. ([2]). The number 
of operations that are needed to build these matrices scales 
as 0{N (^ ). This scaling makes clear that the QML method 
can treat only a limited number of pixels. Therefore in the 
context of an all sky observations it can be applied only at 
modest resolution. 



2.2 Fiducial spectra 

For our fiducial model, we assume the concordance ACDM 
model, with parameters derived from the WMAP7 best fit. 
Thoughout this work, we assume an equation of state of the 
dark energy fixed at w — —1. With these assumptions, it 
is straight forward to calculate the expected power spectra 
Cf ^ and Cf « : 



^GG 



dk 



Cr"=4^ / ^^\k)lf 



(fc), 



Cj^ = 47r 



^A^(fc)jr-(fc)/f(fc), 



(7) 



(8) 



respectively. A^(fc) is the logarithmic matter primordial 
power spectrum, and the filters of the galaxy density dis- 
tribution [If) and the ISW (//^^) are given by: 



lY{k)= I dzb{z)^5M{k,z)ji{kx{z)) 



7/^^(fc) = -2 I dze-^^j,{kx{z)) 



(9) 



(10) 



Here, ^ is the redshift distribution of the galaxy survey in 
question, and we have implicitly used the fact that the den- 
sity contrast in the galaxy survey tracks the matter density 
contrast as: 

fc(n) = / dzb{z) —— 5m {n, z) . 



(11) 



It is well known that the late ISW-LSS cross-correlation de- 
pends not only on the matter fiuctuations on large scales, 
but also on how these are related to the observed galaxy dis- 
tribution, determined by the the product b{z)dN/dz. This 
can be simultaneously estimated using the measurement of 
C( , also exploiting the QML method. 



2.3 Numerical Improvements 

For the reasons discussed above, the QML method is quite 
computationally expensive and prohibitive at high resolu- 
tion. We discuss here some changes which can improve the 
numerics and decrease substantially the execution time with 
a negligible loss of accuracy. 

The predicted Cj is generally non-zero, and its mea- 
surement is the primary object of the ISW measurements. 
However, it is expected to be relatively small, even for the 
largest scales, so it is a good approximation to assume 
Cj — for the fiducial model, which is used to build 
the covariance matrix. Further, the noise matrix N may 
be assumed to be uncorrelated between the CMB and the 
galaxy measurements. Under these assumptions, the Fisher 
matrix becomes block diagonal and the three spectra Cf^, 
Cj , Cf can be estimated independently from each other. 
This reduces the computation cost of the Fisher matrix 
by 50% with respect to the problem with the full covari- 
ance. Moreover estimating just Cj the computational cost 
of the problem decr e ases b y a further factor of 1/6, as in 
IPadmanabhan et all (|2005l '). 

In order to apply the algebra of the QML method, de- 
scribed in Eqs. (|H3|) . one must build the covariance matrix 
C in pixel space and the Fisher matrix F in ^ space. The 
latter is the most expensive task at computationally, largely 
because it requires the inversion of the pixel space covari- 
ance matrix C. This inversion can also introduce numerical 
errors since its eigenvalues naively span several orders of 
magnitude. 

To bypass this issue, we have used inversion-routines 
only on numerically homogeneous blocks thanks to the fol- 
lowing expressions. Given a general matrix A in block form. 



A 



All Ai2 
A21 A22 



(12) 



where An and A22 are non-singular square matrices, then 
it can be shown that the inverse of A is 



A-' = 
with 



-Bll -B11A12A22 1 Qg\ 

-A22 A21BII A22 + A22 A21BIIA12A22 



Bii = (All - ^12^22^^42 



(14) 



For our purposes, we partition the TT, TG and GG 
blocks of C, so that An is the covariance related to the 
CMB temperature sector and yl22 relates to the covariance 
of the galaxy sector. Thus, assuming a fiducial model with- 
out any cross-covariance simplifies the inversion calculation 
significantly. This technique is also applied to the Fisher 
matrix inversion in multipole space (with An — F]jF), ob- 
taining a much better precision with respect to the brute 
force inversion. 
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3 VALIDATION WITH SIMULATED MAPS 

In order to test our implementation of the QML method, 
we created simulated galaxy count maps and CMB tem- 
perature anisotropies follow i ng th e recipe described 



iBoughn. Crittenden fc Turok (119981) fse e also lBarreiro et alj 
l|2008l ') and iGiannantonio et all (120081)'). We e mploy the 



HEALPix Fl program svnfast (|Gorski etahl 120051 ) which al- 
lows one to create airn such that 



(« 



Y' 



) 



C'l Sifi 5„ 



(15) 



where Y, Y' — T, G. The total map for the CMB anisotropies 



is simulated as the sum of three different maps 



T 



ISWc I ISWu I prim 



(16) 



where a^'^ represents the fully correlated ISW effect with 
the galaxy distribution, a^^" is the uncorrelated part of the 
ISW effect and a^™ is the primordial CMB signal. These 
amplitudes are given by 



ISWc _ ^ 



c, 



TG 



yep' 



,ISWu 



= Cm Cfw 



(CJG)2 



(^ISW 



prim _ g /TytT 

In addition for the galaxy count maps we consider 



-iay/C 



GG 



(17) 



(18) 



(19) 



(20) 



where ^'s are Gaussianly distributed complex random num- 
bers, with zero mean and unit variance. They are the seeds 
of the simulations and satisfy {^a^^') ~ ^aa' ■ In this way it 
can be shown that 



; T T »\ _ ^TT 



GG 



; T G*\ _ f^TG 



(21) 
(22) 
(23) 



We have tested the QML approach using these Monte 
Carlo simulations. In particular, we have performed 3000 
realizations for CMB and LSS correlated maps at the 
HEALPix resolution of A''sido = 32 0. For the multipoles, 
we consider the range Ai = [2,95]; i.e., up to the Nyquist 
freque ncy 3A^Bide — 1. Th e standard ACDM cosmological 
model (jLarson et al.ll201ll ) is assumed, as well a survey char- 
acter istics similar to the NVSS catalogue (jCondon et al.l 
119981 ). namely: a similar sky coverage (see next Section), 
a g alaxy den s ity nu mber distribution per redshift given by 
the lHoet all (|2008l ) model, and a bias b = 1.98. 

These simulated maps show that our QML implementa- 
tion leads to unbiased and minimum variance results when 
considering the realistic case of a masked sky, as can be 
seen by comparing the simulations to the projected errors 
from the Fisher matrix. Importantly, we confirm that the 
method is unbiased and minimum variance when the signal 
covariance matrix is block diagonal, i.e. when fiducial cross 
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Figure 1. Tlic average estimates for the Monte Carlo validation: 
the upper and middle panels show the TT and GG auto-spectra, 
respectively, and the lower one is the TG cross-spectrum. We 
compare results for three cases: using realistically masked maps 
without noise in the LSS maps (thick error bars), using full sky 
maps with NVSS-like shot noise (solid line error bars), and as- 
suming both masked maps and NVSS-like shot noise (light dark 
error bars) . We can see that average power spectra from the QML 
all agree very well with the underlying fiducial theoretical power 
spectra (blue lines) . The error bars change according to the noise 
level in the LSS map and the fraction of the sky considered. The 
dark lines are the average of the anafast estimates, which are 
slightly biased at high £ in the two auto-spectra. 



^ http://healpix.jpl.nasa.gov/ 

^ The number of pixels A'^pix is related to the parameter N^^^^ 

through TVpix = 12Af2.^^. 
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Figure 4. APS estimates from WMAP 7 year and NVSS data for the 2.5mJ flux cut. The upper panels show the auto-spectra, TT and 
GG, the lower arc the cross-spectra, TG, in the no binning case (left panel) and in the binning case (right panel). The blue solid lines 
are always the fiducial power spectra. The TG cross-spectrum in the lower right panel is binned with A£ = 9 to make clearer the QML 
estimates respect to the fiducial power spectrum (blue solid line) and the CT^ = (black solid line). We can note the QML estimates 
for the GG auto-spectrum arc slightly larger than the expected fiducial power spectrum. 



power spectrum Cj is set to zero: with the latter approx- 
imation, no difference can be appreciated by eyes on the 
QML estimates and a very small difference can be seen in 
the likelihood constructed by Fisher, which will shown for 
our application to real maps of WMAP 7 yr and NVSS in 
Sect. 

It is important to notice that, while on these large-scales 
the noise contribution in WMAP and future (Planck) CMB 
temperature maps is so low that the CMB noise N can be 
safely neglected, this is not necessary true for large scale 
structure surveys. Depending on the number of sources used 
as large scale tracers, the galaxy density map could be sig- 
nificantly affected by Poissonian noise, which must be taken 
into account. 

The results from the Monte Carlo validation are sum- 
marized in Fig. [l] the upper, middle and lower panels show 
respectively the average estimates for the TT, GG and TG 
spectra derived from the Monte Carlo simulations. Three 
different scenarios are considered, all of which provide un- 
biased averaged estimates in good agreement with the fidu- 
cial model (blue lines), and they differ only in their error 
bars. The first case corresponds to a masked sky (account- 
ing for the NVSS sky coverage and the WMAP KQ75 mask 
with negligible Poissonian shot-noise contribution to the LSS 
map (given by the thick error bars); second, a full-sky case 
with a shot-noise like the that expected in NVSS (see next 



Section for more details) when only sources above 2.5 mjy 
are taken into account (solid line error bars); and, finally 
a more realistic situation where both, the incomplete sky 
and the shot-noise are included in the analysis (light dark 
error bars). The error bars increase when the noise level in 
the LSS map rises and when the fraction of the sky consid- 
ered is reduced, the latter falling approximatively with the 
\//sky/s^y. as expected. 

For comparison, the plots also include (dark lines) the 
average anafast estimation for the full-sky case (dark lines), 
based on the simple HEALPix FFT tool. As it can be seen, 
the anafast estimation is slightly biased at high i in the two 
auto-spectra. 



4 APPLICATION TO WMAP 7 YEAR AND 
NRAO VLA SKY SURVEY DATA 

In this Section we describe the application of our QML code 
to estimating the cross-correlation spectrum between the 
WMAP 7-year CMB maps and the NRAO VLA Sky Survey 
(NVSS) data. 
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Figure 2. WMAP 7 yr (top) and NVSS (bottom) maps at 
HEALPIX resolution N^^^^ = 32 used in this analysis, with the 
respective masks. In the displayed NVSS map, the threshold flux 
is 2.5mjy and the corrections for systematics in declination has 
been applied. 



Figure 3. Comparison of the source redshift distribution used in 
Ho ct al. (2008) (blue) and the one with a redshift dependent bias 
with assuming the CENSORS source distribution (red). 
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Figure 5. Comparison of the binned GG and TG spectra for 
three thresholds in flux, i.e. 2.5 mj (stars with thick error bars), 
5 mJ (diamonds with grey error bars) and 10 mJ (cross with 
error bars). The results appear largely consistent, while the errors 
increase as the number of sources decrease. 



4.1 The maps 

For WMAP data we make use of publicly available prod- 
uctqj- In particular, clean maps at the V and W frequency 
bands have been co-added, using a weighting procedure that 
accounts for the instrumental noise variance per pixel. These 
frequency map s have been clean ed following a template fit- 
ting approach ([Gold et al.ll2011^ . and are those used by the 
WMAP tern to perform c osmological tests, suc h as con- 
straining non-Gaussianity (|Komatsu et al.M201ll ). The co- 
added map has been degraded from its original A^side = 1024 
down to NsidB = 32, since the angular scales associated to 
this resolution (« 2°) is enough to capture almost all the 
signal in the CMB-LSS cross-correlation expected from the 
ISW effect. Following this, the WMAP KQ75 Galactic mask 
(similarly degraded) is applied to the co-added map, in order 
to mitigate the unavoidable foreground contamination in re- 
gions within and near the Galactic plane, and also to remove 
known and intense extragalactic objects such as the Magel- 
lenic clouds and large clusters near the northern Galactic 



^ http://lambda.gsfc.nasa.gov/ 
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Figure 6. Comparison of tiie binned GG and TG spectra wlien 
tlie declination correction is included (diamonds) or not included 
(cross), assuming a flux cut of 10m Jy. 



pole. Finally, the remaining monopole and dipole moments 
outside the mask have been estimated and remo ved. 

The NVSS catalogue jCondon et all 1 19981 ') is a radio 
sample at 1.4 GHz produced with the Very Large Array. It 
covers ~ 85% of the sky, up to an Equatorial declination 
of 5 > —40°. The original survey accounts for ~ 2 x 10® 
sources with fluxes > 2.5mJy. This survey has been widely 
use d in the context of th e ISW studies. It was first used 
by iBoughn fc CrittendenI (|2002| ) to probe the CMB-LSS 
cross-correlation with the COBE data, and a few years 
afterwards it was successfully used by the same authors 
with WMA P data, in the first work re porting such cross- 
correlatio n (IBoughn fc Crittendenll2004l '): this was soon fol- 
lowed by (jNolta at al. II2004I ') with a similar analysis by the 
WMAP team. 

The survey has a somewhat inhomogenou s sensitivity as 
a fun ction of the equatorial declination (see ICondon et al.l 
119981 . for the details), resulting in the mean galaxy density 
that artificially varies with the declination. Therefore, some 
pre-processing is needed in order to mitigate this large-scale 
effect. One of the procedures used in the literature consists in 
defining iso-latitude bands (in equatorial coordinates) and 
imposing that these bands to have the same mean galaxy 
density. In our case, this pre-processing consists of selecting 
first the sources above a particular fiux cut, and then defin- 



Figure 7. Comparison of the binned GG and TG spectra for the 
fiducial model (black) and the one with a redshift dependent bias 
with the CENSORS source distribution (red). Note how the esti- 
mates are stable with respect to the change in the fiducial model, 
and how an evolving bias could alleviate the tension between the 
theoretical GG predictions and the NVSS estomates. 



ing nine bands of equal area, imposing the same mean galaxy 
density number for each band. Finally, we rotate to Galac- 
tic coordinates to compare to WMAP, and then pixelise to 

a HEALPix resolution of N^ide = 32. 

Previous works (e.g., iNolta at al. I |2004| : 



20061 ) 



have shown 



IVielva. Martinez-Gonzalez fc Tuccil 
that the particularities of this pre-processing do not affect 
significantly the results and we confirm this below. We also 
repeat the analysis for different thresholds in flux, namely 
2.5 , 5 , lOmJy, as the higher flux thresholds should be less 
sensitive to possible declination systematics. 



4.2 Source redshift distribution 

To interpret the results of our measurements, we must as- 
sume some redshift distribution dN/dz and potentially red- 
shift dependent bias h{z) for the sample. Given a redshift 
distribution, the average bias can be estimated from the 
measured QML estimates for Cf^; however, here we exploit 
previous measurements for the NVSS sample. 

Historically, the redshift distribution was b a sed o n 
models of the sources by iDunlop fc Peacockl (|l99(J ). 
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and a time-independent b ias of 1.6 was derived by 



iBoughn fc Crittenden (120021'). A larger time- i ndepe ndent 
bias was found by iBlake. Ferreira fc Borrilll (|200J), al- 
beit with a different redshit distribution th an used by 
iBoughn fc CrittendenI (|2002l ) . In lHo et al.l (120081 ) , a new red- 
shift distribution was derived based on a F distribution fit 
which was constrained to give the cross-correlations mea- 
sured between the NVSS survey and SDSS LRG subsam- 
ples: 



dN^ 



dz 



^a + \ 



F(a) 



(24) 



where z. = 0.79 and a = 1.18 0. IHo et all (|2008l ') also esti- 
mates an effective, redshift independent value for the bias as 
b[z) = 1.98. Finally, we also explor e the most recent gal axy 
redshift distribution proposed by Ide Zotti et al.l (|201(J ): a 
fourth order polyno mial fit to the CENSORS distribution 
(|Brookes et al.ll2008l ): 

— ; — = 1.29 + 32. 37^ - 32.89.2^ + 11.13z^ - 1.25/ . (25) 
dz 

A comparison of the two redshift distributions based on Eqs. 
(|24p and ()25|l is shown in Fig. ^. For the second redshift 
distribution in Eg. ([25ll we consider a redshift de pendent 



bias (JMatarrese et al.lll997l : lMoscardini et al.l[l99l ): 



b{z) = bo + 



bi 



D^{z) 



(26) 



where D( z) is the linear gro wth factor in a ACDM Universe. 
Following IXia et all (|2010al ). we choose bo = 1.1, 6i = 0.6, 
7=1. Below we focus on the latter two distributions, and 
examine how the uncertainties impact the derived cosmo- 
logical constraints. 



4.3 Measurements of the spectra 

In Fig. |4] we present the TT, GG and TG spectra obtained 
by our QML up to ^ = 64 (= 2Nside) for the 2.5 mj flux cut 
in NVSS data. Since the signal-to-noise for unbinned TG 
power spectrum is rather poor, we present also the binned 
power spectrum C^'^ over A£ = 9. The binned estimates are 
simply the average of the unbinned estimates inside the bin. 
For plotting purposes, we associate ■>/) ]„ 



for the uncertainty in the binned estimate. Unless otherwise 
stated, all the maps have been corrected for the declination 
effect. 

To investigate potential systematic problems, we com- 
pare the dependence of the TG and GG spectra on the dif- 
ferent threshold fluxes for NVSS considered here in Fig. [5] 
(TT is not shown since it is of course unchanged.) Overall 
the APS estimates agree very well when varying the flux 
threshold, with larger error bars for larger flux threshold, as 
expected given the fewer objects and resulting larger Poisson 
errors. In Fig. |S] we examine the importance of the correc- 
tion for declination systematics in N VSS for a flux cut o f 
10 mJ. Our resuh for GG agrees with lBlake fc Walll (|2002D . 
confirming that with a conservative flux cut of 10 mJ the 
declination systematics in NVSS is negligible. 

* Note that we have corrected t he norma l izatio n factor 
of the r distribution assumed in lHo et alj (120080 ■ 



Finally, in Fig. [12] we compare the GG and TG results 
obtained by assuming the two redshift distributions in Eqs. 
()24l) and (|25p , including also the dependence of the bias on 
the redshift in the latter case. The QML estimates are sta- 
ble with respect to such different physical assumptions in the 
two fiducial power spectra. By considering the more physi- 
cal scenario in which the bias evolves in redshift, the tension 
between the NVSS auto-spectrum and the theoretical pre- 
dictions could be alleviated. 

We note few important findings in our estimates of the 
angular power spectrum of the WMAP 7 year - NVSS cross- 
correlation. First, the estimate of the GG power spectrum 
in the NVSS map is slightly larger than our fiducial model. 
Our esti mates for the NVSS auto-p o wer sp ectrum agree very 
well with lBlake. Ferreira fc Borrilll (|2004l ). who used an op- 
timal estimator similar to ours on a NVSS map of the same 
resolution of the one used here. The stability of the Ci 
estimates with respect t o diffe rent flux threshold found in 
iBlake. Ferreira fc Borrilll (|2004l ) is also very similar to what 
we find. 

As yet it is unclear whether this deviation could be 
caused by some systematics in the NVSS data or should 
be ascribed to a genuine physical effect, as a n effective 
bias la rger than ~ 2, which is usually assumed. IXia et al.l 
(|2010aj) estimated a larger discrepancy at lower multipoles 
and explained this effect as result of non-negligible pri- 
mordial non-Gaussianity, caused by the large-scale scale- 
dependence of the non-Gaussian halo bias. However, the 
value inferred for the coupling non-Gaussian parameter 
/jvL is much l arger than the limit s imposed by CMB 
analyses (e.g., iKomatsu et al.l 1201 ll : ICurto et al.l 120111 ). 
The /jvi c onstraints derive d from the CMB-LSS cross- 
correlation (|Xia et al.ll2010a ) provide lower values, in better 
agreement with the CMB tests. In addition, these authors 
also showed that when other LSS data sets are used ( in par- 
ticular, the QSOs sample of the SDSS lRichardsll2009l ') . such 
non-Gaussian deviation is not found. 

We also note an estimate lower than the fiducial model 
in the first bin of the TG spectrum, less than 2a (as obtained 
by the Fisher matrix) lower than the fiducial model. This 
low value was n ot obtained in the previous investigation by 
iHo et all (|2008D . 



5 DARK ENERGY CONSTRAINTS 

In this Section, we constrain the dark energy density Qa 
using the information contained in the ISW-LSS cross- 
correlation power spectrum, estimated through our QML. 
We assume the errors on the measured Cj are Gaussian, 
and calculate the relative likelihoods of SIa using 



-21n[£(fiA)] = x'(nA) 
where 



2 
Xmiii ■ 



(27) 



x'i^A) 



[^- 



TG, obs 

e 



Cri^A)] C-,} [C,T^'°'^= - Cj^'^iQA)] . (28) 



Here C^ '° ° are the unbinned estimates of the cross- 
correlation power spectrum, and Cf [Qa) are the theoreti- 
cal predicted power spectrum. The matrix Cui is the covari- 
ance matrix between different multipoles, which allows for 
correlations among non-diagonal terms which arise in the 
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Figure 8. The likelihood for Ha obtained by the Fisher prescrip- 
tion in Eq. JgOj, with the 95 % and 68 % C.L. for the threshold 
flux of 2.5 mj (blue), 5 mj (red), 10 mj (black) in NVSS, respec- 
tively. 



Figure 10. Comparison of likelihood contours for Ha. obtained 
by the Fisher prescription in Eq. II29I I when accounting (dashed 
line) and when not accounting (solid line) for shot noise in NVSS 
data. The threshold flux in NVSS has been chosen as 2.5 mJ. 




Figure 9. Comparison of likelihood contours for Qa obtained by 
the Fisher prescription (black solid line), the covariance computed 
by the Monte Carlo for CMB only (red short dashed line) and for 
both CMB and LSS (blue dashed line). The threshold flux in 
NVSS has been chosen to be 2.5 mJ. 




Figure 11. Comparison of likelihood contours for Ha obtained 
by the Fisher prescription in Eq. ||29|I when considering the full 
covariance (red, solid line) and when using the approximation of a 
block diagonal signal covariance and Fisher matrix. The threshold 
flux in NVSS has been chosen as 2.5 mJ. 



Amin 



is the minimum value of x with 



presence of masks, 
respect to 0,^- 

We compare the hkehhoods obtained by different pre- 
scriptions for the covariance matrix. The first prescription is 
to use the unbinned QML estimates and the Fisher matrix 
as its covariance matrix: 



C[e, = {F- 



(29) 



An alternative prescription is to construct the covariance 
matrix C by averaging over Monte Carlo realisations of the 
maps. For every model Ha, we can define the covariance C 
with A'' simulated CMB and LSS maps 



N 



E 



[cn(«A) 



c/^(nA)][c/,^(nA)-c/,^(nA)] 

iV 



(30) 
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Figure 12. Comparison of likelihood contours for Qa obtained by 
the Fisher prescription in Eq. I I29I I for the two choices of redshift 
distributions: solid for Eq. II24I I and dashed for Eq. I I25II . The 
threshold flux in NVSS has been chosen as 2.5 nij. 



where the Cj ^ are the estimates for every single realization 
i and the Cf is their theoretical value. However, we expect 
the covariance matrix not to depend strongly on the cosmo- 
logical model, one can just consider the case with SIa = 0, 
and since Cj {Q.a = 0) = 0, the covariance becomes. 



pMC _ 



E 



N ■ 



(31) 



We can build Cw in Eq. (|3ip either by using random re- 
alisations of only the CMB maps and the single, true NVSS 
map, or by creating a realisations of both CMB and LSS 
maps. In both cases, w e generate our covariances on 1000 
realisations, as done bv lVielva. Martinez-Gonzalez fc Tuccil 
(J2Q06l l. We also examine how the probability contours for 
JIa depend on the various assumptions such as the thresh- 
old flux cut used for the NVSS map or the sources redshift 
distribution. 

We evaluate the likelihood with the various different 
prescriptions by sampling the x^ on values of ^Ia , < S^a < 
0.95. The other cosmological para meters are k ept fixed to 
the values determined by WMAP (JLarson et a l. 2011) for 
the standard ACDM model. As default NVSS description, 
the Eq. (|24p model is assumed, with a bias of 1.98. 

By adopting the Fisher matrix prescription in Eq. (1291) . 



as tightest constraint we obtain JIa ~ 0-69_„'22 rn 501 ^^ 
l(2)cr confidence level (CL) for the lowest flux threshold of 
2.5 mj, see the blue dashed line in Fig. [51 An Einstein-de 
Sitter Universe is disfavoured at more than 2 a CL for the 
lowest flux threshold in NVSS, consistent with earlier mea- 
surements. Note that the conditional probabilities for SIa 
agree for the different flux thresholds considered. 

By building the covariance through realizations of the 
CMB maps while keeping the NVSS map fixed, we obtain 
the probability distribution given by the red dashed line of 
Fig. [91 We find ^a = O-G^tl'llU'll] ^t l(2)cr CL. Using in- 



stead the covariance derived from realizations of both CMB 
and LSS maps, the probability distribution given by the blue 
dashed line of Fig. Ewe find Oa = O-^Sto^^l'io'AJ) at 1(2)0" 
CL. Note that the constraint based on the Fisher covariance 
is tighter than the one based on a Montecarlo covariance 
keeping fixed the NVSS map, but looser than the Monte- 
carlo covariance obtained with CMB and LSS uncorrelated 
maps. Overall, the three likelihood prescriptions are consis- 
tent, although some of the differences might be ascribed to 
the (unexplained) discrepancy between the C^'^ estimates 
and the theoretical predictions, on which the Montecarlos 
are based. 

Given the agreement among the three different likeli- 
hood prescriptions, we can use the Fisher prescription for the 
covariance to test other dependences of the analysis. In Fig. 
(|10[) we verify the importance of taking into account the shot 
noise in the NVSS map: by not taking into account the shot 
noise the probability contours for Oa would be much tighter, 
even for the maps with the most sources. We then study the 
impact of approximating the signal covariance matrix (and 
consequently the Fisher matrix) as block diagonal, i.e. con- 
sidering Cf — for the fiducial underlying model. This ap- 
proximation is not essential for our app roach, whereas it is 
necessary for lPadmanabhan et al.l (|2005|) in which only Cf 
is estimated. As mentioned in Sect. 2.3, the difference in the 
power spectrum estimates is barely visible and we have ver- 
ified by MonteCarlo that this approximation does not alter 
the optimality of the method. On the real data considered 
here. Fig. (|lip shows how the constraints with the full covari- 
ance and Fisher are a bit tighter than thos e with the block 
diagon al assumption. As already noticed in lCruppuso et al.l 
(|2009h . conditional probability slices are much more sensi- 
tive to small changes than QML estimates. 

In Fig. p2p we compare t he redshift distributio n esti- 
mated with CENSORS data bv lde Zotti et all (|2010l ) in Eq. 
(|25|) . including a red shift dependent bias as from Eq. (), with 
the one adopted bv lHo et al.l (|2008l ). with an effective bias 
of 1.98. 



6 DISCUSSIONS AND CONCLUSIONS 

We have developed an optimal estimator for the angular 
power spectrum of the cross-correlation between CMB and 
maps of large scale structure, which in parallel estimates 
their auto-spectra. This has been tested using an ensemble 
of randomly generated maps, and we have demonstrated the 
robustness of the QML estimates for the TT, TG and GG 
power spectra. Our QML implementation extend similar op- 
timal estimato rs limited only to the galaxy auto power spec- 
trum (Blake, Ferreira fc B orrill 2004 ) or only to the c ross- 
correlation power spectrum (jPadmanabhan et al.ll2005l ). 

We have applied our method to WMAP 7 year and 
NVSS data, the best public data sets at present for study- 
ing the ISW cross-correlations. Our method makes no as- 
sumptions, and allows to measure the cross-correlation with 
optimal errors and to exploit the full cosmological informa- 
tion contained in the maps, though our analysis is limited 
to a pixel resolution off. 8°. While the NVSS map contains 
known declination systematics, we correct for these and find, 
as has earlier work, that they appear to have little effect on 
the measured cross-correlations. In agreement with previous 
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studies, we detect a non-zero cross-correlation, and have also 
seen a slight excess in the NVSS auto-angular power spec- 
trum compared to that expected theoretically. 

We have translated these measurements into the quanti- 
tative constraints on the fraction of dark energy in a ACDM 
model which can be obtained only by the cross-correlation 
of WMAP and NVSS, estimating Qa while keeping fixed all 
the other cosmological param eters to the WMAP 7 yr best- 
fit values ([Larson et al.ll20ffl ). We have compared three dif- 
ferent prescriptions for estimating the covariances: using the 
Fisher matrix computed by our QML, on Monte Carlo relisa- 
tions of the CMB maps only and creating Monte Carlo real- 
isations of both CMB and LSS maps. We have found a good 
agreement among the Oa probability contours obtained from 
these three different likelihood prescriptions. The width of 
this probability contour depends mainly on the flux thresh- 
old and associated level of Poisson noise in the NVSS map, 
but the signal amplitude seems largely independent of the 
flux. The constraint from the likelihood prescription based 
on the Fisher matrix we derive from the cross-correlation be- 
tween WMAP 7 yr and NVSS data is Qa = OmtlU'l^fJ) 
at 1(2)0" confidence level (CL) for the lowest flux threshold 
of 2.5 mj. Such value is quite consistent with the concor- 
dance cosmology This result agrees with that expected from 
a typical survey with sk y fraction and noise property as the 
NVSS, and agrees with lVielva. Martinez-Gonzalez fc Tuccil 
l|200a ), but is somewhat w eaker than the one obtained by 
the no n-optimal analysis bv lPietrobon. Balbi and Marinuccil 
l|2006l ) based on needlets. It is not clear if this discrepancy 
is due to the lower resolution considered here or the ne- 
glection of shot noise in the NVSS m ap in the analysis by 
IPietrobon. Balbi and Marinuccil (|2006l ). 
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